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We consider the "modified Minimal Analytic" (mMA) coupling that involves an infrared cut to 
the standard MA coupling. The mMA coupling is a Stieltjes function and, as a consequence, the 
paradiagonal Pade approximants converge to the coupling in the entire Q 2 -plane except on the 
time- like semiaxis below the cut. The equivalence between the narrow width approximation of the 
"q^" ■ discontinuity function of the coupling, on the one hand, and this Pade (rational) approximation 

' of the coupling, on the other hand, is shown. We approximate the analytic analogs of the higher 

f"*) , powers of mMA coupling by rational functions in such a way that the singularity region is respected 

■ by the approximants. Several comparisons, for real and complex arguments Q 2 , between the exact 
| i ' and approximate expressions are made and the speed of convergence is discussed. Motivated by the 

, success of these approximants, an improvement of the mMA coupling is suggested, and possible uses 

■ in the reproduction of experimental data are discussed. 

' PACS numbers: 

I. INTRODUCTION 

In perturbative QCD (pQCD) calculations, the coupling a(Q 2 ) = a s (Q 2 )/n shows non-physical singularities at low 
energy (— q 2 = Q 2 < 0.1 GeV 2 ). The aim of analytic QCD (anQCD) is to give a coupling Ai(Q 2 ) which is analytic 
at low Q 2 (> 0) and reproduces the high energy behavior of a(Q 2 ). Using the Cauchy theorem we can write QCD 
' running coupling in the integral form 

8: °(o*)=if a) 

where r\ — > +0 and p\{p) is the discontinuity function of a(Q 2 ) along the cut axis in the complex Q 2 -plane at n-loop 

■ approximation given by p^ 1 i '\o~) — Im[a'™~^'(— a — ie)]. 

Q\ ' The Minimal Analytic (MA) procedure of Shirkov and Solovtsov [l[ removes the pQCD contribution of the unphysical 
. cut, < — a < A 2 , keeping the discontinuity elsewhere unchanged 

■ A (MA) (0 2 ) - - r da (2) 

h ; 

This expression doesn't have singularities for Q 2 > and, as required, reproduces the high energy behavior of a(Q 2 ). 

In fact, A[ MA \Q 2 ) is analytic in the entire complex plane of Q 2 with the exception of negative semiaxis, reflecting 
the analyticity properties of the space-like observables. 

We will consider a modification to the MA coupling [l[ of Eq. Following Ref. 0, the lower limit of integration 
is increased to a certain value Mq which we set to be ~ M 2 (~ 10 _2 -10 _1 GeV 2 ). Thus, the modified MA coupling 
is given by the following dispersion relation: 

4 mMi) (Q 2 ) = -f^. (3) 

One of the motivations for modifying the MA in this way is to include the point Q 2 = and its vicinity in the 
analyticity region of Ai(Q 2 ), something reflected by analyticity properties of the space- like observables, among them 
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the V-type Adler function, where the infrared (IR) cutoff Mq is: Mg ~ M 2 (~ A^j-,). 1 When M5 2 > 0, the analytic 
coupling A\ is a Stieltjes function, with the radius R of convergence for its Taylor expansion around Q 2 = given by 
R = Mq. Stieltjes functions have the attractive property that their paradiagonal Pade approximants must converge 
to such functions at any complex value of the argument (except on the cut) when the order index of the approximants 
increases. Such a behavior is not guaranteed in MA model because there Mq = and, consequently, Q 2 = is not a 
point of analyticity and A^ is not a Stieltjes function. 

In this work we investigate the behavior of the coupling _4( mMj4 ) (Q 2 ) (with M 2 ~ M 2 ) of the mMA model ([3]) and 
its paradiagonal Pade approximants. In Sec. [ill the basic definitions of rational (Pade) approximants and of Stieltjes 
functions are presented. In Sec. IIIII we demonstrate that the narrow width approximations for the discontinuity 
function pi(er) of A^ M are equivalent to approximating Ai M (Q 2 ) by paradiagonal Pade approximants; we 
investigate the convergence of the latter approximants to the "exact" function values A^™ MA \Q 2 ) for Q 2 > 0. In 
Sec. IIVI we extend our analysis to the higher power analogs An MA ^ of a n (n > 2), construct the approximants for 
j^mMA) k age( j on the aforementioned paradiagonal Pade approximants for j^ lM ^ anc [ investigate their behavior 
for Q 2 > 0. In Sec. [V] we investigate the behavior of the aforementioned approximants for complex values of Q 2 . In 
Sec. I VII we propose an extension (improvement) of the mMA model, motivated by the success of the narrow width 
approximations of the discontinuity function p\{cr) in mMA; further, we point out the prospects for applications of 
such simple analytic QCD models to fitting experimental data. Section [VIII represents the summary of our results. 



II. PADE APPROXIMANTS AND STIELTJES FUNCTIONS 



Consider a function f(z), of a complex variable z, with power expansion about the origin given by 

N+M 

f(z) « ]T f n z n . (4) 

71=0 

The Pade (or: rational) approximant, R^ t (z), is defined as the rational function 

r n M {z) = , (5) 

satisfying the condition that its expansion about z = matches N + M + 1 terms in the power series expansion in 

M 



dH) . The paradiagonal Pade R M 1 can be written (using partial fractions) in the following way: 



<- 1 w = E^ w 

& -|- 4> n 

n—1 

where z n are zeros of the denominator in Eq. ([5]). 

If the function f(z) is a Stieltjes function with a finite radius of convergence R > 0, i.e. 

m =ir dT ^L, (7) 

7T J l+TZ 

with g(r) a nonnegative function, thus f(z) having the real coefficients of its series expansion around z = given by 

(-IY r 1/R 

fn = / r n g(r)dT, (8) 



then a strong convergence theorem in the theory of Pade approximants applies which says: in the limit M — > oo, 
Rm +J (z) (J > —1) is equal to f(z), and the poles of R^ +J (z) are simple poles which lie on the interval of the 
negative real axis given by — oo < z < —R (Ref. Q, Sections 5.2, 5.4). 

In particular, when f(z) is a Stieltjes function, the Pade i?^ _1 (z) of f(z) has the form ([6]) with d n > and z n > 0, 
the poles —z n having the following ordering: — oo < —zm < —zm—i < ■ • • < —z\ < —R. 



We note that Mo is a renormalization scheme (RSch) dependent quantity, just like A is (e.g., A in V-RSch and MS RSch differ); 
so = A/q/A 2 is RSch-independent. Therefore, while we expect Mq ~ M n; Mq so) is a free parameter of the model. 
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III. APPROXIMATION OF THE COUPLING BY RATIONAL FUNCTIONS 

In order to establish a relationship between A^ nMA \Q 2 ) and a Pade approximant, we first note that the modified 
MA coupling in Eq. §5§ is a Stieltjes function, as defined in Eq. (O by identifying: Q 2 — z, a = 1/r and pi(cr) = rg(r). 
It has the convergence radius M5 2 , and the series expansion about Q 2 = given by 

4— > w >, = ±f M V^ = f;L„(<m (9) 



?i=0 



L n = [ —^~ / r"3(r)dr, (10) 



where 



and g(r) = ap\(a). For simplicity, we want the expansion coefficients to be dimensionless. To fulfill this, we 
introduce dimensionless variables s — er/A 2 = 1/t and u = Q 2 /A 2 in Eq. ©, where A 2 = Aq CD (in MS convention; 

Aq CD ~ 10 _1 GeV 2 ). Consequently, Eqs. © and (10]) obtain the form 

A[ mMA \uA 2 ) = i / dtM- = J2 L,y\ (11) 



with 

r (-1)' 



1/so 

* n fl (t)dt, (12) 

where u = Q 2 / 'A 2 = z/A 2 , s = Mq/A 2 , L„ = Z„(A 2 )" are dimensionless, and g(t) = p 1 (A 2 /t)/t. 

An important consequence of the fact that j^™- MA ^ j s a Stieltjes function is that the series in Eq. @ only converges 
within the finite disc in the Q 2 -pl&ne (Q 2 < Mq), elsewhere the series is divergent. However, we will see that, 
because of the aforementioned theorem Q , only few coefficients of this divergent series are needed in order to evaluate 
A < ] mMA \Q 2 ) in regions well beyond the convergence disc, via approximate analytic continuation in the form of Pade 
approximants R 1 ^ 1 ~ 1 (z). 

Following Ref. Q, we can approximate 9(s) = g(t = l/s) = spi(a = A 2 s) by a singular series expansion 

SP1 ( S A 2 ) _ g(l/s) _ 6(s) W 



— «E/^( S " S «), (13) 



7r 

n— 1 



where s = a/A 2 , s n = M 2 /A 2 and f 2 are all positive dimensionless quantities. In Ref. [J], this type of approximation 
was applied to the spectral function of the vector channel vacuum polarization function iIy(Q 2 ). It was motivated 
there by the highly singular behavior of the vector spectral function in the limit N c — > oo. In our case we aren't 
working in that limit, but we note that the spectral function pi(er) of A^-[ nMA \Q 2 ) is positive; although Q(s) has 
finite values, such a function can be well approximated as a sum of positively weighted Dirac deltas for the purpose 
of integration. We will see that this approximation is equivalent to approximating the mMA coupling A < ] mMA \Q 2 ) 
as the Pade approximant i?^ _1 (u). 

Using the singular approximation (|13p for 9(s) in Eq. (jlip . we obtain 

n=\ n=l ^ 11 

where F 2 — f 2 A 2 . This relation is the same as R^~ 1 (z) Pade of j(™- MA ) given in Eq. (J6j) , with z = u, z n = s n and 

d - f 2 

n— 1 n— 1 " 
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Therefore, by approximating 9(s) by positive Dirac deltas we obtain ^4( mMj4 ) written in a (paradiagonal) Pade form. 
Further, the following relation holds [3j: 



\A[ mMA) (uA')-RZ~ L (u)\ < K 



M-l, 



^/s + u- y/s^ 



+U+ y/sS 



2 M 



(16) 



where if is a constant. Therefore R 1 $ [ ~ 1 {u) — > _4( mMj4 ) (uA 2 ) when M — + oo. 

In our evaluations we will use the perturbative coupling given by the solution of the two-loop renormalization group 
(RG) equation 



i^ L \uk 2 ) = 



-1 



m(i + w T i(z ± )y 

where W^i(z) is the Lambert function (branches n = ^l), and the argument z± is given by 



z± 



o i(±7r-0/3 o /ci) 



(17) 



(18) 



where Q 2 = \Q 2 \e^>, \u\ = |Q 2 |/A 2 , (3 = (11 - 2n//3)/4, ci = /VA), ft = (102 - 38n / /3)/16. For < (f> < n, 
the branch W-i(z+) is chosen, for — n < <j> < 0, the branch W+i(z_) is chosen. At low energies (\Q 2 \ < 10 GeV 2 ), 
the number of active quark flavors is nt — 3. For details see Ref. []|. With the coupling given in Eq. (fT7| . the 
discontinuity function becomes 



^• ) ( S A 2 )=Im 



-1 



ci(l + W+i(z_(s))) 



(19) 



with = — 7r and |u| = a/ A 2 = s in Eq. (|T8|) . These expressions are in fact valid at any n-loop level (n > 2) in the 't 
Hooft renormalization scheme. 
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FIG. 1: A[ MA) {Q 2 ) (continuous), and A{ mMA) (Q 2 ) at two-loop level for various values of Mo cut. Throughout this paper we 
use rif — 3 and A^jg- = 0.35 GeV. 

We set A = 0.35 GeV (at n/=3). 2 This value can be changed later when we fit experimental data. Here we will 
compare numerically the accuracy of the rational approximants with the "exact" numbers, i.e., those obtained by 
evaluating numerically integrals Eq. ([3]). 



2 In Ref. @|, the value A (n/=3) ~ 0.4 GeV was obtained in MA by requiring that the MA model reproduce measured values of QCD 
observables at higher energies (|Q| > 10 GeV) in MS scheme. 
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TABLE I: The first twenty coefficients computed with the discontinuity function at two-ioop given in Eq. I|19[> . for Mo = 2M 7V . 
In parentheses, the results for Mq = 3M,r are given. All the numbers here are presented as rounded off at the sixth digit. 







L n 




Lo 


0.151418 (0.127410) 


Lio 


0.270477 (6.95492 • 10" a ) 


Li 


-0.0408795 (-0.0149290) 


Ln 


-0.384736 (-4.40027 • 10~ 5 ) 


Li 


0.0353764 (0.00582799) 


Ll2 


0.551687 (2.80615 • 10" 5 ) 


L 3 


-0.0380741 (-0.00281121) 


Ln 


-0.796468 (-1.80155 • 10~ 5 ) 


U 


0.0453111 (0.00149464) 
-0.0571406 (-8.40628 ■ 10~ 4 ) 


1/14 


1.15653 (1.16323 ■ 10" 5 ) 


U 


Lis 


-1.68780 (-7.54792 ■ 10" 6 ) 


u 


0.0748244 (4.90454 ■ 10" 4 ) 


Lie 


2.47387 (4.91884 ■ 10" 6 ) 


L 7 


-0.100605 (-2.93630 ■ 10~ 4 ) 


Ln 


-3.63999 (-3.21771 ■ 10" 6 ) 


L 8 


0.137943 (1.79193 ■ 10~ 4 ) 


Lis 


5.37405 (2.11200 ■ 10~ 6 ) 


U 


-0.192008 (-1.10983 ■ 10~ 4 ) 


-1/19 


-7.95836 (-1.39043 • 10~ 6 ) 



In Fig. Q] we compare mMA coupling (for various values of the IR cutoff Mo) with the standard MA coupling 3 at 
low positive Q 2 's. These values are calculated by performing integrals in the dispersive relations ([3]) and ([2]) for each 

Q 2 . An effect of increasing the value of the cutoff Mq in Eq. ([3]) is to decrease the values of A < f nMA \Q 2 ) at low Q 2 . 

We can compute L n coefficients of mMA appearing in Eq. (fTTj) by using the discontinuity function Q19p in integrals 
(fT2)l . All the integrations are done numerically. The first twenty coefficients are shown in Table |U With these 
coefficients we can compute the paradiagonal Pade approximants up to Rf (u) (for higher Pade's we need 

more coefficients). 

All the rational approximants i?^ _1 (u) up to Rf Q (u) for A^C lMA \uK 2 ) are shown in the Figs. [21 for two choices 
of the IR cutoff in Eq. ([3]): Mo — 2M T , 3M^. We recall that Mq is the convergence radius of the Stieltjes function 
A^ nMA \Q 2 ). Since this radius is larger in Fig. [2Kb) , the convergence of the Pade approximants is also faster there. 

The rational function i?f (Q 2 ) matches A ( ™ MA) {Q 2 ) in Fig.UJa) quite well for all < Q 2 < 10 GeV 2 , and in Fig.^b) 
the agreement is even better. 

In order to give an idea of how the approximants work, we list in Table [TT] the values of Q 2 at which the deviation 
from Ai 1lMA \Q 2 ) reaches 0.1% (the deviation increases with increasing Q 2 ), for the case Mq = 2M n . For the 
Pade approximant with M — 10 the deviation is less than 0.1% for all Q 2 < 2.5 GeV 2 , and less than 

1% for all Q 2 < 5.4 GeV 2 (when the Pade index is M = 20, these values increase to 10.33 GeV 2 and 22.63 GeV 2 , 
respectively). We recall that these values are much higher than the convergence radius of the Taylor series, Eqs. © 
and (HB): R = Q 2 = M 2 = 0.0784 GeV 2 . We see that the first pole (-si) converges from below to the IR cut 



3 A\ '(0) = l/Po ~ 0.44; however, the slope dA\ >(Q 2 )/dQ 2 at Q 2 = is infinite. 
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TABLE II: The values of Q 2 = uA 2 at which the Pade approximants R\\ 1 (up to R\ ), for M = 2M V , reach 0.1% deviation 
from ^4^ mMA '. Further, the first few values s n = M 2 /A 2 and f 2 are shown. 



r>M — 1/ \ 
% W 


Q 2 [GeV 2 ] 


Si 


it 


S2 


/I 


S3 


fi 


R { ! 


0.01 


3.704 


0.56085 










Rh 


0.08 


0.912073 


0.0263116 


13.2502 


1.62407 






R.i 


0.20 


0.747835 


0.00947914 


1.62259 


r\ nr on An r 

0.0538095 


30.2372 


3.19243 


R-i 


0.38 


0.698643 


0.00498205 


1.04773 


0.0185788 


2.62744 


0.0831117 




0.60 


0.677037 


0.00309638 


0.871413 


0.00968867 


1.46119 


0.0275631 


Re 


0.88 


0.665567 


0.00211868 


0.7913 


0.00604702 


1.11052 


0.0140762 


Rl 


1.20 


0.658731 


0.00154383 


0.747386 


0.00417231 


0.951402 


0.00869082 


Ra 


1.58 


0.654321 


0.00117625 


0.720466 


0.00306913 


0.863709 


0.00596662 


R<> 


2.02 


0.651309 


0.00092658 


0.702677 


0.00236037 


0.809515 


0.00438199 


Rio 


2.50 


0.649158 


0.000749085 


0.690269 


0.00187583 


0.773381 


0.0033713 



— s — — Mq/A 2 = — 4Af 2 /A 2 (—so — —0.64 in our case) when index M increases, as expected. For example, when 
M = 10, we obtain: si = 0.6492 (ra s = 0.6400), s 2 = 0.6907, etc. When the cut M is increased from 2M V to 3AU 
(=> -s = -1.44), we obtain for M = 10: s x = 1.4605 (» s = 1.44), s 2 = 1.5525, etc. 



IV. APPROXIMANTS FOR THE HIGHER POWER ANALOGS OF THE ANALYTIC COUPLING 

To obtain the analogs of the higher powers of the analytic coupling in this formalism, we use relations given in Ref. 
At the 3-loop level truncated series for the analytic coupling, we have: 

A 2 (Q 2 ) = A 2 (Q 2 )- Cl A 3 (Q 2 ), A 3 (Q 2 ) = MQ 2 ), (20) 

where 

An[Q ' ffi-\n-l)\ d(lnQ2)«-i- W 

The correspondence between the powers a k of the perturbative coupling a(Q 2 ) = a s (Q 2 )/n and the above quantities 
is: a k i— > Ak- The couplings Ak are the analytic versions ("analogs") of higher powers needed for evaluation of 
observables. We note that in general Ak ^ A k (for further discussion, c.f. Sec. Ill of Ref. @)- 

Using the dispersive integral expression Eq. ([3|) in Eqs. pi"]) , we obtain explicit expressions for Ak's in terms of 
integrals of the (perturbative) discontinuity function p\ 



A { 2 mMA \uA 2 ) = JL r ^Q- ds - ClA i mMA \uA% (22) 

fa* j so {s + uy 



We see from here that the quantities A^ MA \Q 2 ) have the same location of the singularities as A^ 71 ^ 1 ^ : Q 2 < —M'q 
(= — A 2 so). On the other hand, the poles Q 2 = —M 2 = — s n A 2 (n = 1, . . . , M) of the Pade approximants R^ 1 for 
Ai, Eq. (|15p . do reflect, in a discretized manner, the singularity cut of A^ 1 ^ 1 ^] namely, these poles appear on the 
negative real axis in the range: — oo < — M| f < —M 2 I _ l < • • • < — Mf (< Mq). Such a pole structure of i?^ _1 (w) 
is guaranteed because j\ L ^ nMA ^ is a Stieltjes function Q. On the other hand, the analytic higher power analogs 
are no longer Stieltjes functions. Therefore, if we construct the Pade approximants 1 (u) 

for A^ MA \Q 2 ) (fc > 2) in the usual way, i.e., based on the first 2M coefficients of their Taylor expansion around 
Q 2 = 0, there is no longer the guarantee that all the poles of the obtained rational functios lie on the negative real 
axis (and below —Mq). For instance, the Pade approximant Rf computed from the series of A3 has a non-physical 
pole on the positive real axis. One way to avoid this problem is to compute the rational approximants of life's directly 
from the derivatives of the Pade approximant of ^4^ mMj4 ' ) . If i?j^ -1 (u) are the paradiagonal Pades that approximate 
A^ nMA \uA 2 ), we shall define the rational approximants R 2 {u) and Rs(u) for A 2 mMA \uA 2 ) and A 3 lnMA \uA 2 ), 
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respectively, in the following way: 



where 



i? 2 (u) = B,2(u) — cxRs(u), R 3 (u) = R 3 (u), 
R n (u) 



( _ 1} n-l d n-l R M-l (u) 



PS~ l (n-i)\ a(inu)"- 1 



(24) 



(25) 



A comparison between Rk and their respective Au, computed from R\ Q {u) (i.e., Pade index M — 10), is shown in 



0.005 



0.004 



o.os: 



0.000 




0.0003 



-0.0002 



-0.0004 



-0.0006 




2 4 6 

O 2 [GeV 2 ] Q 2 [GeV J ] 

FIG. 3: Comparison between: (a) A^™ MA ^ (continuous) and f?2 (dashed and dot-dashed), for two choices of IR cut Mo = 2M W 
and Mo = 3M n , respectively; (b) the same as in (a), but now for .A3 M and R-j,. Pade index is M = 10. 

Figs. [3] The accuracy of these approximants decreases notably compared to the R^~ X approximants of A < ( nMA ^ . The 
accuracy of the approximants increases when the cut Ma increases, e.g., from 2M W to 3Af w . 

V. COMPLEX ARGUMENTS 



In evaluation of observables, sometimes evaluation of the analytic coupling and their power analogs at complex 
values of arguments is needed (e.g., see App. C of Ref. Q). For some complex arguments Q 2 = \Q 2 \ exp(i(f>) [u = 
Q 2 /A 2 = \u\ exp(i(f>)], the results based on the approximant Rf (u) (Pade index M — 10), are shown in Figs. |4][7] 
When (j) — it/4, Figs. H] (a), (b) show that the approximant for A^™ MA \Q 2 ) works with less than 1% error up to 
\Q 2 \ ~ 10 GeV for the real part, and \Q 2 \ ~ 3 GeV for the imaginary part. The approximants R2(u) and i?3(it) for 

A2 and A3, (see Figs. EH?]) fail already at lower \Q 2 \. Figs. 0(a), (b) show the approximant for A^" lMA \Q 2 ) when 
4> = 'in /A. Figs. HE indicate that the relative accuracy of our approximants decreases: (I) when the index n of the 
power analog A n increases; (II) when the considered ray Q 2 = \Q 2 \ exp(i0) comes closer to the time-like semiaxis (\(f>\ 



closer to tt). Aspect (I) can be understood from the fact that, by our construction, R n involves 



R 



M-l 
M 



1) derivatives 



of the Pade approximant R\ 
singularity cut of AnS. 

The accuracy can be increased if we compute R\ 

. (mMA) 



Aspect (II) is also to be expected, because rays with 4> ~ ±tt are close to the 



~ Rm 1 as a higher order Pade (higher index M). For instance, 



Figs.[5](a), (b) show the approximants for with cf) — 7r/4, which are calculated from i?2o( u ) (i- e -j Pade index 

M = 20). The latter is calculated from the first 40 coefficients of the Taylor series of4 mMA) K 2 ). We see that in 
the case of Mo = 2M W , the deviation of Re[i?2] from ReL42] becomes discernible to the eye (1% deviation) only for 
Q 2 > 5 GeV 2 ; this is to be compared with Fig. [6] (a) (where Pade index M = 10). If we used = instead (see 
Figs. [3] where M = 10), the deviations of Re[R 2 ] from Re[*4 2 ], for M = 20 and M = 2M 7T , on the scale of Figs. M 
would turn out to be discernible to the eye only for Q 2 > 6 GeV 2 , i.e., the convergence turns out to be even better 
than in the case <fi = n/4. The deviations of Im[i?2] from Im[^2] in Fig. [5] (b) are not discernible to the eye. We also 
see that the deviations of R2 from A2 cannot be seen by the eye in Figs. [5] when the cut value Mq increases to SM^- 4 



At the level M = 20, it is important that the 40 Taylor coefficients L n and the coefficients of the Pade approximant 



be calculated 
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Both aforementioned aspects (I, II) that decrease accuracy are, however, not very important in practice when 
evaluating observables. Namely, the higher order contributions (~ A n ) are very suppressed in anQCD (even when 
\Q 2 \ is low), and the contributions of Q 2 near the time-like axis in the contour-type of integrations (e.g., for the 
semihadronic r decay ratio ry) are usually suppressed by the integrand. 





VI. PROSPECTS OF APPLICATIONS IN FITTING EXPERIMENTAL DATA 

In Eq. © we considered the dispersive relation for A\ with an IR cr-cutoff cr cu t = Mq (~ M 2 ) imposed on the 
perturbative QCD discontonuity function pi(c) = Im[a(Q 2 = — o — ie)}. Such a cutoff results in the analyticity of 
Ai around Q 2 = 0, thus reflecting the analyticity of space-like observables V(Q 2 ) in the complex plane excluding the 
time-like semiaxis but including a regime around Q 2 = 0. 



to high accuracy (to at least about 20 and 30 digits, respectively), in order to avoid numerical instabilities connected with cancellation 
of large numbers. 
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FIG. 7: (a) Real parts of .A3 M and 7?3 with complex arguments Q 2 = \Q 2 | exp(i7r/4); (b) same as in (a), but for imaginary 
parts. 

In MA model, the scale A can be fixed so that it reproduces measured values of QCD observables at higher energies 
Q > 10 GeV A„ /= 3 sa 0.4 GeV) j|6f. However, then MA gives too low value of the massless strangeless (AS = 0) 
semihadronic T-decay ratio: ry ~ 0.14 The experimentally measured value of this quantity is: r T = 0.204±0.005 

Q. The latter value can be reproduced in MA with A(„ /=3 ) w 0.4 GeV only if the current masses of light quarks 
(m u ,md ~ 1 MeV) are replaced by much larger (constituent) masses (m Ul md ~ 0.25 GeV) [10( and the threshold 
effects become very important. 

By introducing IR cutoff er cu t = Mq , the coupling A\ gets further diminished at low Q 2 , and thus further diminishes 
the value of r T . To remedy this situation, we can, in the simplest way, simulate the unknown behavior of pi(tr) (Fig. [9]) 
at a < Mg by adding a simple positive Dirac delta peak: 5pi(sA 2 ) — 7r/^ 1 <5(s — s_i), where < s_i < s a (= ^/q/A 2 ). 
This would then allow us to achieve, in the model, the correct value of r T while still maintaining the analyticity of 
A\ around Q 2 = 0. Thus, the full discontinity function in such a "delta-modified" MA model (dmMA) is 

^ dmMA) (,sA 2 ) = 6( S - So )MsA 2 ) + tt/V( S - S _i) , (26) 

where Q(x) is the Heaviside step function (+1 for x > 0, zero otherwise), and pi(sA 2 ) is the pertubative QCD 
discontinuity function: pi(sA 2 ) = a(—sA 2 — ie). This leads to the following A\\ 

A[ dmMA \uK 2 ) = I Hds Pl{sA2) + / " 1 . (27) 
7r7 Sn S + U U + S_i 
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Re(A 2 ) 
M = 20 



3M» 




-0.0010 



iei 2 [Gev-j 

FIG. 8: Now 7?2 is computed from the Pade R20 of A^"" 1 ^ (M = 20) : (a) real parts of A 2 mM A ' (continuous) and R2 (dashed 
and dot-dashed), for the IR cut value Mo = 2M^, 3M^-, respectively, with complex arguments Q 2 = \Q 2 \ exp(i7r/4); (b) same 
as in (a), but for imaginary parts. For Mq = 3M„-, no deviations of R2 from A^™ 1 ^ 1 ^ can be seen by the eye. 



Applying the Pade 1 (it) approximation to this analytic coupling we obtain 

1 ("A ) W i? M („)= ^ _= ^ ^ 

n=-l(n/0) n=-l(n^Q) 

This has the same form as the Pade applied to _4^ mMA ' ) (MA with IR cut), Eq. (fT5)) . but just with one more 

term (n = —1). This model has three dimensionless model parameters: so,s_i,/_i. All are positive and ~ 1. As 
presented, the model is considered in the 't Hooft scheme (J3j — for all j > 2). The scale parameter A n/= 3 is fixed 
by fitting the model to experimental values of observables at high energies (|Q| > 10 GeV), such as T decay, e + e~ — > 
hadrons, Z — » hadrons. The values of low energy QCD observables (\Q\ ~ 1 GeV), such as r r and Bjorken polarized 
sum rule, are sensitive to the values of parameters sq, s_i and f—i; therefore, the latter are to be fixed by fitting to 
the experimental values of such observables. 




FIG. 9: This Figure shows the IR cut so = Mq/A 2 (when Mo = 2M^) for perturbative pi that was used starting in Eq. Q. The 
unknown low momentum part (0 < s < so) can be simulated by adding a Dirac delta, say irf'L 1 8(s — s_i), where < s_i < sq. 



11 



We see that even with this modification (dmMA) of the MA, the evaluation of the couplings A n (uA 2 ) (n — 1, 2, . . .) 
is made simple and efficient by using the Pade approximant (|28[) ; for the evaluation, it suffices to know the (three) 
parameters sq, s_i and /_i and the first few coefficients Lk [Eqs. (fTT1l - (fT2|) and TableHj. 



VII. SUMMARY 



We worked with the minimal analytic (MA) model modified (mMA) by an IR cutoff <r cu t = Mq ~ M 2 for the 
perturbative discontinuity function. In such a model, the analytic coupling Ai(Q 2 ) is analytic in the entire Q 2 - 
complex plane excluding the semi- axis Q 2 < — M 2 . The analytic properties of such (mMA) coupling reflect the 
analytic properties of space-like QCD observables T>(Q 2 ), among them analyticity in the point Q 2 = and its vicinity. 
Further, such a (mMA) coupling is a Stieltjes function of Q 2 . This implies that it will be efficiently approximated by 
(paradiagonal) Pade approximants R 1 ^ I ~ 1 (Q 2 ), i.e., R^~ 1 (Q 2 ) converges to Ai(Q 2 ) at any point of analyticity when 
index M increases (M -> oo). The coupling A X {Q 2 ) in the form of R%~ 1 (Q 2 ) can be easily and efficienty evaluated 
[i.e., without performing time-consuming dispersion-type integrations © for each Q 2 ] just by knowing the first few 
coefficients of Taylor expansion of Ai(Q 2 ) in powers of Q 2 . 

We showed that for real and complex arguments the paradiagonal Pade approximants of the analytic coupling 
Ai MA \Q 2 ) are precise at low positive Q 2 's. This high precision range of positive Q 2 's increases fast when the order 
index M of the Pade approximant increases. When Q 2 's are complex, the precision range of |Q 2 |'s decreases when 
Q 2 approaches the singularity cut. The analytic analogs An MA ^ of higher powers a n = (a s /vr) n (n > 2) are then 
evaluated as combinations of logarithmic derivatives of the approximant R^ I ~ 1 (Q 2 ) [w A\(Q 2 )]. The approximants 

obtained in this way for A^" M show less precision when n increases and/or when Q 2 approaches the singularity cut. 
These approximants still work fine if we increase the order index M of the Pade. However, high precision is needed 
only for the n = 1 case, because the higher couplings A^h nMA ' 1 get strongly suppressed (even at low \Q 2 \) in analytic 
QCD when n increases. Further, when evaluation of observable involves contour integration (such as, for example, in 
the case of r T ), the contributions of A n (Q 2 ) get supressed by the rest of the integrand when Q 2 comes close to the 
singularity cut. While we generally used for the IR cutoff a cut (~ M 2 ) the specific value o- cut (= -M5 2 ) = 4M 2 , we also 
showed that the conclusions in this work are independent of the specific value chosen, by comparing various results 
for M = 2Mn and M = 3M». 

We further suggested an inclusion of one additional Dirac delta function to the mMA discontinuity function 9(a — 
M 2 )pi(a) at low energies where the precise behavior of pi is unknown - Dirac modified MA model (dmMA). Such a 
modification maintains the analyticity at Q 2 — and its vicinity, and allows us to reproduce the experimental value 
of the semihadronic r decay ratio r T . Such a modification keeps the same form of the Pade approximants R^^ 1 (Q 2 ) 
of j\(d mMA ) as j n m MA case and allows us to evaluate them (and the higher power analogs) in an easy and 
efficient manner. The (three) parameters of such a model can be determined by requiring that the model reproduces 
the measured values of low energy QCD observables [ill ]. 
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